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We present a new, fully analytical point scattering model which can be applied to arbitrary anisotropic 
magneto-electric dipole scatterers, including split ring resonators (SRRs), chiral and anisotropic plasmonic scat- 
terers. We have taken proper account of reciprocity and radiation damping for electric and magnetic scatterers 
with any general polarizability tensor. Specifically, we show how reciprocity and energy balance puts constraints 
on the electrodynamic responses arbitrary scatterers can have to light. Our theory sheds new light on the mag- 
nitude of cross sections for scattering and extinction, and for instance on the emergence of structural chirality in 
the optical response of geometrically non-chiral scatterers like SRRs. We apply the model to SRRs and discuss 
how to extract individual components of the polarizability matrix and extinction cross sections. Finally, we 
show that our model describes well the extinction of stereo-dimers of split rings, while providing new insights 
in the underlying coupling mechanisms. 



I. INTRODUCTION 

Research in the field of metamaterials is driven by the pos- 
sibility to control the properties of light on the nanoscale by 
using coupled resonant nanoscatterers to create optical ma- 
terials with very unusual effective medium parameters. En- 
gineering arbitrary values for the effective permittivity e and 
permeability fj, would allow new forms of light control based 
on achieving negative index materials 1-3 , or transformation 
optics media 4 that arbitrarily reroute light through space. In 
order to reach such control over e and /i, many metamaterial 
building nanoblocks have previously been identified as hav- 
ing an electric and magnetic response to incident light, in- 
cluding split ring resonators (SRRs) 5-10 , rod-pairs 11 , cut-wire 
pairs 12 , fishnet structures 13-15 and coaxial waveguides^. In 
many instances, the nanoscatterers are not only interesting as 
building blocks in subwavelength lattices of designed e and 
fi. The building blocks are in fact very strong scatterers with 
large cross sections 17-20 , comparable to the large cross sec- 
tions of plasmonic structures. Therefore, metamaterial build- 
ing blocks are excellently suited to construct magnetic anten- 
nas, array waveguides and gratings in which electric and mag- 
netic dipoles couple and form cooperative excitations, in anal- 
ogy to the functionality imparted by plasmon hybridization 21 . 
Experiments outside the domain of effective media have ap- 
peared only recently. These experiments include experiments 
by Husnik et al^, and Banzer et alr^ that quantify the extinc- 
tion cross section of single split rings under differently polar- 
ized illumination, as well as a suite of experiments on coupled 
systems. These experiments include extinction measurements 
on split ring dimers 23 that point at resonance hybridization, as 
well as reports of magnetization waves 24 , structural and geo- 
metrical chirality in arrays, as evident in e.g. massive circular 
dichroism 25-32 , and chiral effects in split ring stereo-dimers 
studied by Liu et a/. 33 . 

In order to understand the light-metamatter interaction in 
systems of strongly coupled magneto-electric scatterers, it is 
important to understand how individual metamaterial building 
blocks are excited and how they scatter. So far, explanations 
of the observed phenomena have mainly rested on two pil- 



lars. On the one hand, data are compared to brute force finite- 
difference time-domain (FDTD) simulations of Maxwell's 
equations, usually showing good correspondenc o 5 i 6 i 8 i 17 i 1 8 i 25 i 26 . 
While the FDTD method is a rigorous method, such numeri- 
cal experiments do not in themselves provide insight into how 
split rings scatter or hybridize in coupled systems. There is 
general consensus that to lowest order, metamaterial inter- 
actions in lattices of scatterers like SRRs must be described 
by magneto-electric point-dipole interactions. Hence, simple 
LC circuit models with dipolar coupling terms are the second 
main interpretative tool to predict, e.g., frequency shifts due 
to electric and magnetic dipole-dipole interactions in lattices 
and oligomer a 9 ' 23,33 ' 34 . To rationalize this LC circuit intuition, 
several authors have analyzed current distributions obtained 
by FDTD simulations in order to retrieve the microscopic pa- 
rameters {i.e., the polarizability) underlying such a dipolar in- 
teraction model, and in order to estimate multipolar correc- 

While there is general consensus that to lowest order, meta- 
material interactions must essentially be magneto-electric 
point dipole interactions, we note that the dipolar circuit mod- 
els in use so far have some significant shortcomings. Strictly 
speaking, the electric circuit theories lack the velocity of light 
c as a parameter. Hence, they contain no retardation or in- 
terference, they violate the optical theorem, do not predict 
quantitative cross sections, and fail to predict the effects of 
super- and subradiant damping on resonance linewidths. A 
fair comparison of intuitive point-dipole ideas with actual data 
is therefore impossible, unless a fully electrodynamic theory 
for magneto-electric point dipoles is derived. Such a theory 
would generalize the electric point scattering theory that is 
well known as very effective means to describe random media, 
extraordinary transmission and plasmon particle arrays^ 9 -^. 
In this paper we derive exactly such a theory for general 
magneto-electric scatterers. We show how reciprocity and en- 
ergy conservation restrict the full magneto-electric response 
via Onsager constraint a 42 ' 43 , and a new magneto-electric op- 
tical theorem for the full polarizability tensor. This tensor 
not only includes an electric (magnetic) response to electric 
(magnetic) driving, but also off-diagonal coupling in which a 
magnetic (electric) response results from electric (magnetic) 
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driving. While our theory sheds no light on the microscopic 
origin of the polarizabilityii, we show how electrodynamic 
polarizability tensors can be directly constructed from LC cir- 
cuit models. Furthermore we predict how extinction measure- 
ments and measurements of radiation patterns (i.e., differen- 
tial scattering cross section) can be used to quantify the polar- 
izability tensor. 

The paper is structured in the following way: Firstly, in 
Section HI] we derive the general theory, taking into full ac- 
count reciprocity, the optical theorem and radiation damping. 
In Section [Til] we apply this theory to set up the polarizability 
of the archetypical metamaterial building block, a single SRR. 
In Section [TV] we show which set of experiments can be used 
to retrieve the tensor polarizability a. We find that magneto- 
electric coupling directly implies circular dichroism in the 
extinction of single split rings, evidencing the utility of our 
theory to describe structural chirality^ 5 -^. Thirdly, we show 
in Section [V] that the theory can be simply applied to obtain 
quantitative scattering spectra of coupled systems. By way of 
example we examine the case of two coupled resonators in the 
stereodimer configuration reported by Liu et alM. 



the incident light 

a =( aEE aEH ) . (2) 

Here, olee quantifies the electric dipole induced by an applied 
electric field. The tensorial nature of olee is well appreci- 
ated in scattering theory for anisotropic particles, such as plas- 
monic ellipsoids 48 . By analogy with the electric response to 
electric driving quantified by olee, the tensor olhh quantifies 
the magnetic dipole induced by a driving magnetic field. Fi- 
nally, the off-diagonal blocks represent magneto-electric cou- 
pling. The lower diagonal olhe quantifies the magnetic dipole 
induced by an incident electric field, and a eh quantifies the 
electric dipole induced by an incident magnetic field. Such 
magneto-electric coupling is well known to occur in the con- 
stitutive tensors of metamaterials 7 . Indeed, the first metama- 
terials consisted of split ring resonators, in which there is a 
magnetic response without any driving magnetic field in nor- 
mal incidence experiments 49 . However, the relative strength 
of magneto-electric coupling in the polarizability, i.e., oleh, 
and olhe have not been experimentally quantified. 



II. MAGNETO-ELECTRIC POINT SCATTERER 



B. Electrodynamic Onsager relation 



A. Polarizability 

A paradigm in scattering theory is the point dipole scat- 
terer— -41,45,46 to model scattering by very small, but strongly 
scattering particles. Generally, incident fields E and H in- 
duce a (complex) current distribution in an arbitrary scatterer. 
It is the express point of this paper to assess what the scat- 
tering properties are of subwavelength scatterers with strong 
electric and magnetic dipole moments, as this represents the 
physics expected of metamaterial building block s 9 ' 23 ' 33 i 34 ' 47 . 
Therefore we retain only electric and magnetic dipole terms, 
neglecting higher order multipoles. In such a theory, each 
scatterer is approximated as an electric dipole with an electric 
dipole moment p = oleeE that is proportional to the driv- 
ing electric field E. The proportionality constant is the po- 
larizability oleE' In this paper, we derive a generalized point 
scattering theory for metamaterials that includes a magnetic 
dipole moment m on an equal footing with the electric dipole 
moment p. In the most general case, the electric dipole mo- 
ment p and magnetic dipole moments m are induced by both 
the external electric and magnetic fields E and H according 
to 
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m 



= a 



(1) 



Throughout this paper we suppress harmonic time depen- 
dence e~ ZUJt . We use a rationalized unit system that signif- 
icantly simplifies all equations and is fully explained in Ap- 
pendix lAl In Eq. (fl]i, a is a 6x6 polarizability tensor, which 
consists of four 3x3 blocks, each of which describes part of 
the dipole response to the electric or magnetic component of 



There are several constraints on a. In addition to any sym- 
metry of the scatterer itself that may impose zeros in the po- 
larizability tensor, these constraints are due to reciprocity and 
to energy conservation. We start by examining the constraints 
imposed by reciprocity. It is well known from the field of 
bi-anisotropic materials 43 that reciprocity imposes so-called 
Onsager constraints on the most general constitutive tensors 
relating (D, B) to (E, H). Already Garcfa-Garcfa et al£L 
proposed that such Onsager constraints carry over directly 
to electrostatic polarizabilities. Here we rigorously derive 
Onsager relations for electrodynamic magneto-electric point 
scatterers. By, definition, the electric and magnetic fields due 
to the induced p and m are equal to 



E on t 
H ou t 



= G°(r,r') 



P 



(3) 



with a dyadic Green tensor G° that describes the field at posi- 
tion v = (x, y, z) due to a dipole at r' — (x' , y', z'). The 6x6 
Green dyadic of free space can be divided in four 3x3 blocks 



G°(r,r') = 



G% H (r,r>) 



G HE {ry) G% H (ry 



(4) 



The 3x3 diagonals correspond to the familiar known elec- 
tric field Green dyadio 4 ^ 4 ^ and magnetic field Green dyadic 
of free space, which in our unit system (see Appendix) both 
equal 

G EE (r, V) = G° HH (r, r<) = (Ik 2 + VV)^^. (5) 

The off diagonal blocks correspond to the mixed dyadics that 
specify the electric field at r due to a magnetic dipole at r', 
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respectively the magnetic field at r due to an electric dipole at 
r'. Explicitly: 



G° EH (r,r') = -G° HE (ry) 

( d : 
= ik I -d z a 

d y -d x 



d z -d y \ ik\r-r'\ 



(6) 



In this work we focus solely on scatterers made from recip- 
rocal constituents, as is commonly true for the metallic scat- 
terers that constitute metamaterials. Since the materials that 
compose our scatterers (typically gold and silver) satisfy reci- 
procity microscopically, the polarizability tensor must also 
lead to a scattering theory that satisfies reciprocity. 

To derive reciprocity constraints on a, it is sufficient to ex- 
amine the Green function in the presence of just one point 
scatterer at the origin. This Green function that quantifies the 
field at r 2 due to a source at ri in presence of a single scatterer 
at r s can be written a o 39 ' 40 ' 50 

G{n,r 2 ) = G°(r u r 2 ) + G°(r 2 , r s )aG°(r s , n), (7) 

Reciprocity requires for any Green function G (similarly split 
in four blocks) that 



G EE {r 2 ,ri) G EH {r 2l ri) 
G H E{r 2 ,ri) G ffff (r 2 ,ri) 

G EE {r 1 ,r 2 ) -G EH (ri,r 2 ) 
-G HE (ri,r 2 ) G H H{ri,r 2 ) 



(8) 



which is equivalent to noting that swapping source and detec- 
tor leaves the detected field unchanged, up to a change in sign. 
Specifically, Lorentz reciprocity requires a transpose for the 
diagonal 3x3 blocks, meaning that swapping a detector and 
source of like character leaves the detected field unchanged. 
An extra minus occurs for the off-diagonal terms, i.e., when 
swapping a magnetic (electric) detector with an electric (mag- 
netic) source. It is easy to verify that Eq. (0 is indeed satisfied 
by the free space Green function G°. 

Using this fact, we evaluate Eq. (0 for the Green function 
in Eq. (0 to find if reciprocity constrains a. Since reciprocity 
is clearly satisfied for the first term in Eq. (Q, we now focus 
on the second term 

G^r^r^aG^r^n) = G {r lirs )aG {r Sl r 2 ). (9) 

Expanding the matrix products in Eq. (0 while making use 
of the reciprocity of the free Green function results in the On- 
sager relations for the dynamic polarizability: 

a EE = a EE , a HH = a r HHl and a EH = -c*# B 

(10) 

These relations are identical in form to the Onsager relations 
for constitutive tensors 43 , but are now derived on very differ- 
ent grounds for the polarizability of electrodynamically con- 
sistent point scatterers. This gratifying result shows that the 



general point dipoles proposed in this work can be used as mi- 
croscopic building blocks for an exact scattering theory that 
describes the formation of bi-anisotropic media from dense 
lattices of scatterers in the effective medium limit. Indeed, 
since the point scattering building blocks fulfill the Onsager 
constraints, they are natural building blocks to derive effec- 
tive media constitutive tensors by homogenization that also 
satisfy the Onsager relations. 



C. Optical theorem 

It is well known in point scattering theory for electric 
dipoles that polarizability tensors are not solely limited by 
reciprocity and spatial symmetry, but also fundamentally by 
energy conservation. Indeed, energy conservation imposes 
an 'optical theorem' that constrains the polarizability of an 
electric dipole scatterer to ensure that (in absence of mate- 
rial absorption) extinction equals scattering 40 . We proceed to 
examine these constraints imposed on a. Let us first recapitu- 
late the well known case of a scalar electric scatterer^i .. An 
electric scatterer will absorb and scatter part of the incom- 
ing light, that together make up the extinction of a dipole. 
Extinction for an electric scatterer corresponds to the work 
done by the incident field E in order to drive the dipole 
p. The work per optical cycle needed to drive p equals 
W RcE • Re^j 3>, where <s;g> indicates cycle aver- 
aging. Evaluating the work per cycle, and dividing it by the 
incident intensity I in = \E\ 2 /2Z (Z the impedance of the 
host medium) leads to cr cx t = W/L m — 4nklma EE . Scat- 
tering corresponds to far field radiation radiated by the dipole 
p. According to Larmor, the cycle-averaged scattered power 



is^lP = 
sections 



^^-\p\ 2 . Hence one obtains the well known cross 
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Coxt = 4:Trklma EE and <r scatt = —k 4 \a EE \ 2 . (11) 

Equating extinction to scattering for nonabsorbing particles to 
impose energy conservation, gives rise to the optical theorem 
for the polarizability 



Ima = -fc 3 |aB£;| 2 



(12) 



This equation for instance shows the well-known fact that a 
real (electrostatic) ao, such as Rayleigh's polarizability a = 
3V(e — l)/(e + 2) for a small sphere of dielectric constant 
e, never satisfies the optical theorem 52 . An electrostatic «o 
can be made to satisfy the optical theorem by adding radiation 
dampin g 40 ! 41 to obtain the dynamic polarizability 
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a 



3 



(13) 



It is easy to verify that the albedo of a scatterer with polariz- 
ability given by Eq. ( fT3l is 



Cscat 
Ccxt 



1 + |fc 3 Imo;o ' 
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confirming that radiation damping indeed transforms any loss- 
less electrostatic polarizability (Imao = 0) into a scatterer 
that satisfies the optical theorem. Also material loss included 
in ao via e evidently leads to a lossy scatterer a < 1, as 
expected. Many alternative derivations of Eq. ([T3T > have ap- 
peared, for instance by making a size parameter expansion of 
dipolar Mie coefficients^. 

Inspired by the case of a simple electric dipole, we now 
generalize the optical theorem and the concept of radiation 
damping to the full 6x6 tensorial polarizability of arbitrary 
magneto-electric scatterers. In this case, the work done per 
unit cycle by the incident field E- m and H m to drive p and m 
is equal to 

W= < RcE m • Re^ + RcH m • Rc^ > (14) 
at at 



which evaluates to 

W = ^Hm 



a 



Eu- 
H-„ 



(15) 



where (•)* indicates complex conjugate. The power per solid 
angle radiated by the induced dipoles in a direction specified 
by a unit vector f is easily found by calculating the far-field 
Poynting vector from Eq. ([3]). The result is composed of three 
terms: 



dP 



dP n dP„ 



dn 



— Rc(p x m) ■ f, 



(16) 



The first term in Eq. (fTBT l represents the scattered radiation 
of just the electric dipole p, which integrates to a total scat- 
tered power given by Larmor's formula. The second term in 
Eq. (fT~6b represents the radiation pattern of just the magnetic 
dipole m, again given by Larmor's formula. Note that both 
terms simply represent the well known sin 2 9 donut shaped 
radiation pattern for p and m. The third term, however, can 
completely change the radiation pattern, as it contains the in- 
terference between the fields of p and m. Hence the relative 
phase between the induced p and m is important for the dif- 
ferential scattering cross section. To obtain the total scattered 
power, one should integrate Eq. (fT~6b over all solid angle. The 
interference term integrates to 0, as is easily seen from the fact 
it is an odd function of f. Therefore, Larmor's formula imme- 
diately generalizes, and the scattering cross section equals: 



just the polarizability tensor a. without reference to any in- 
cident field (-Bin, -Hm)- In order to eliminate (E- ln , Hi n ) we 
make the assumption (verified below for split rings) that a can 
be diagonalized. We call the eigenvectors Vi, and denote the 
eigenvalues, which we will refer to as 'eigenpolarizabilities', 
with Ai. Expanding the incident field at the position of the 
origin in the orthogonal eigenvectors 

* i 

and with otVi = A%Vi and (vi\vj) = Eq. (fl~8T) reduces to 



1=1 



A; 



6 
i=l 



(20) 



with strict equality for lossless scatterers. Since this equation 
must be satisfied for any choice of incident wave (i.e., any 
combination of c,), we find a generalized optical theorem for 
6x6 polarizability tensors that can be expressed in terms of 
the eigenpolarizabilities as 



-k 6 \Ai\ > ImAi 



Mi = 1 . . . 6, 



(21) 



again with strict equality for lossless scatterers. Eq. (EH im- 
plies that the polarizability tensor represents an energy con- 
serving scatterer, if and only if each of its 6 eigenpolarizabil- 
ities are chosen to satisfy the simple scalar optical theorem 
(Eq. ( fT2l ) derived for electric scatterers. This general opti- 
cal theorem highlights the importance of two new quantities: 
the eigenpolarizabilities, and the corresponding eigenvectors 
of the point scatterer polarizability. 

In Eq. ( [Oi l we reviewed the well-known addition of radi- 
ation damping required to make electrostatic polarizabilities 
satisfy the optical theorem. Since metamaterial scatterers are 
frequently treated via electrostatic circuit models, it would be 
extremely fruitful to generalize this method to general 6x6 
electrostatic polarizability tensors. It is now evident, that we 
can simply apply the scalar recipe to each eigenpolarizability 
separately. An alternative notation for this method is: 



P 



47T . 

3Z 1 



P 
m 



(17) 



Equating extinction to scattering results in a condition that 
must be satisfied for any incident field (E- ln , H- m ) 



Im 



a 



3 fc 



Ei n H; 



ol a. 



Ei n 



(18) 



Due to the tensorial character of a it is not immediately ev- 
ident how to extract a useful optical theorem that constrains 



-tHl 



(22) 



We note that this expression, which is identical to Eq. ( fT~3b 
upon replacement of l/( ) by matrix inversion, provides a 
unique relation to translate a magneto-/electrostatic polariz- 
ability tensor ao derived from RLC circuit theory, to the cor- 
responding electrodynamic polarizability that satisfies the op- 
tical theorem. We can hence consistently assess how intuitive 
ideas based on a microscopic RLC circuit model for electro- 
static dipoles lead to quantitative predictions for extinction, 
scattering, as well as resonance hybridization, diffraction and 
super/sub radiant damping in coupled systems, such as peri- 
odic systems, or arbitrary finite clusters. 
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FIG. 1. Split ring radiation patterns corresponding to the polarizability tensor eigenvectors. Panel (a): (Sketch) A single split ring resonator 
can have an electric dipole moment p along the x-axis due to charging of the split. Circulating current j in the ring gives rise to a magnetic 
dipole moment m in the z-direction. Panels (b,c): Radiation patterns of the two eigenmodes of an SRR in the case of no off-diagonal magneto- 
electric coupling (r/E = 0.7, r]H = 0.3, r\c = 0). The electric dipole moment oriented along the a>axis radiates most of its amplitude in the 
ky, kz plane, while the magnetic dipole oriented along the z-axis radiates mostly into the kx, ky plane. Panels (d,e): radiation patterns of the 
eigenvectors with magneto-electric cross coupling (r/c = 0.4). Panel (f): indication of the polarization of the light radiated by the eigenvector 
with largest eigenvalue (panel (c)). Light is linearly polarized for wave vectors along the cartesian axes, but elliptically polarized in general. 
The direction of strongest circular dichroism in extinction and scattering is in the xz-plane. 



III. POLARIZABILITY OF SPLIT RING RESONATORS 



A. Symmetry 



The symmetry constraints that set which elements of cxsrr 
are zero, are valid both for the electrodynamic and electro- 
static polarizability of split rings. 



As an example of our general theory we consider the spe- 
cific example of split ring resonators. The electrostatic po- 
larizability of split ring resonators was discussed for instance 
by Garcia-Garcia et al 47 . We consider the LC resonance of 
an infinitely thin split ring in the xy plane, with split oriented 
along the x axis, as shown in Fig.[T}a). Incident electric field 
polarized along the x direction gives rise to an electric dipole 
p = (a E x E E x ,0, 0) oriented along the split of the SRR. As in 
an LC circuit, the charge separation generated over the capac- 
itive split relaxes as a circulating current, hence giving rise to 
a magnetic dipole m = (0, 0, a h f E E x ) in the z direction, in 
response to a driving E-field along x 54 . The same is valid vice 
versa: an applied magnetic field along z induces a magnetic 
dipole moment m = (0, 0, aff H H z ) along the z direction. 
The associated current accumulates at the gap, giving rise to 
an electric dipole moment p = (a x f E H z , 0, 0) driven by H z . 
If we assume that the LC resonance really only involves p x 
and m z , we find that the polarizability tensor is filled with ze- 
ros, except for the four contributions described above. Hence 



B. Quasi-electrostatic RLC model 

We will now construct the electrodynamic polarizability by 
starting from an electrostatic polarizability derived from a sin- 
gle resonant RLC equation of motion. Therefore we take a 
common resonant frequency dependence out of the tensor el- 
ements, writing 



^static 
a SRR 



a{uj) 



( 



\ 



>1E 










-i-qc 



ir\ c \ 






o m ) 



(24) 



where t]e, f]c an d t\h are constant and a(ui) is a Lorentzian 
prefactor 



V 



or 



(25) 







&SRR. = 



ZX 

HE 



*eh\ 



af H 



Here, ujq is the SRR resonance frequency wq 



/LC 



, 7 is the 



damping rate due to the Ohmic loss of gold and V is the phys- 
ical particle volume. As in the plasmonic case, this approxi- 

(23) 

mation is coined 'quasi-static', as it does contain frequency oj, 
but does not contain the velocity of light c. The polarizability 
obtained from the quasi-static polarizability once the radiation 
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damping term is added (section Ull CI Eq. (l22l i) is called 'dy- 
namic polarizability'. In this formulation, all the frequency 
dependence, and the units of cx.srr are contained in a(u>). 
The parameters ije, r\H and r\c are dimensionless. For a loss- 
less split ring t}e,t]h and rjc are all real. We assume that all 
losses are introduced via 7. To determine the sign of t]e, t\h 
and rjc, we expect that for very slow driving the charge (cur- 
rent) on the capacitor directly follows the driving E (H )-field, 
implying tje > and r/n > 0. The sign of rjc follows simi- 
lar reasoning, After charge build-up, charge associated with a 
p x = Ke(a(u!)i]Ee~ lult E x ) relaxes as counter-clockwise cur- 
rent, giving rise to a negative m z = ~Rc(a(uj)irjce~ lut E x ), 
implying that sign 77c = sigmjE- 

C. Limit on magneto-electric coupling 

Having constructed an electrostatic polarizability in accor- 
dance with RLC circuit models proposed in earlier reports, 
we apply radiation damping according to Eq. (l22l to obtain a 
scatterer that has a correct energy balance: 

<*srr = KfAT 1 - \kH. (26) 

So far we have not explicitly discussed absorption loss, except 
through the inclusion of the material damping constant 7 in 
the quasi-static polarizability. Starting from a quasi-static po- 
larizability with quasi-static eigenpolarizabilities Af atlc , the 
albedo for each eigenillumination Vi can be expressed as 



1 1 + |fc 3 Im^f atic ' 

It follows that for any lossy scatterer the imaginary part of 
each eigenvalue ^4| tatlc of the electrostatic polarizability ten- 
sor must be positive to ensure < a < 1. In the case of a 
tensorial ot with loss included as in Eq. (12411 . d25l ). one needs 
to explicitly verify that each eigenvalue has positive imaginary 
part. The eigenvalues of Eq. d24"l i are^4^ atlc = a(w)A± with 

A ± = -. Since lm(a(uj)) > and 

A± are real, we find that both eigenvalues have positive imag- 
inary part only if both A + and A_ are positive. Thus, loss sets 
an additional constraint on the polarizability tensor, and limits 
the magneto-electric coupling to 

|»jc| < Vmm- (28) 

This result implies a very important limitation on magneto- 
electric scatterers: it states that a magneto-electric cross cou- 
pling (r]c) can only be generated if there is a sufficiently 
strong directly electric, and directly magnetic response. We 
note that this constraint is very similar to the constraint on 
the magneto-electric cross coupling in constitutive tensors de- 
rived for homogeneous bi-anisotropic media in Ref. |43| that 
recently attracted attention in the framework of proposals for 
repulsive Casimir forces 55 - 56 . While our derivation was spe- 
cific for split rings, we note that similar constraints hold for 
all magneto-electric scatterers. In the presence of material 



loss, the magneto-electric coupling terms are limited by the 
fact that all electrostatic eigenpolarizabilities must have posi- 
tive imaginary part. 



IV. PREDICTED SCATTERING PROPERTIES OF SINGLE 
SPLIT RINGS 

In the remainder of the paper we discuss some insights that 
the proposed magneto-electric point scattering theory pro- 
vides in how split rings scatter. In this section we will con- 
sider the eigenmodes and the radiation patterns of a single 
SRR for a. given by Eq. d26b . Next, we predict which set 
of experiments will provide full information on the elements 
of the polarizability tensor. We will show how the extinction 
cross sections can be translated back to retrieve SRR polariz- 
abilities and magneto-electric cross polarizabilities of a single 
SRR. Although the results we present are general, we use a 
specific set of parameters for all the figures presented in this 
paper. These parameters are chosen to fit to the properties of 
split rings that are resonant at A = 1.5/zm(o;o/27r = 200 THz 
and that consist of 200 by 200 nm gold split rings with a 
thickness of 30 nm and a gap width of 90 nm. Thus we take 
V = 200 x 200 x 30 nm 3 . We set the damping rate to be 
that of gold 7 = 1.25 • 10 14 s _1 as fitted to optical constants 
tabulated in in Ref. H3- We use tje — 0.7, tjh = 0.3 and 
rjc — 0.4. These parameters were chosen because (A) they re- 
produce quantitatively the extinction cross section under nor- 
mal incidence along the z-axis measured by Husnik et a/. 17 , 
and (B) they fit well to our transmission data on arrays of dif- 
ferent densities of split rings taken at normal incidence 9 and as 
a function of incidence angle 58 . The chosen values correspond 
to on-resonance polarizabilities olee — 4.6V, olhh — 2.1V 
and cleh = 2.5V, all well in excess of the physical SRR vol- 
ume V as is typical for strong scatterers. Finally, we note that 
the calculated albedo fits well to the albedo a — 0.5 to 0.75 
calculated by FDTD by Husnik et al^i. 



A. Radiation patterns and eigenvectors of the polarizability 
tensor 

In Fig. Q] we consider the eigenstates of the split ring po- 
larizability tensor presented in Eq. ( l26| >. We first assume that 
the cross coupling terms are absent, i.e., rjc — 0, in which 
case the polarizability tensor is diagonal, with eigenpolariz- 
abilities a(u))r)E and a(u))rjH. The corresponding orthogonal 
eigenmodes are {p x ,m z ) = (1,0) and (p x ,m z ) = (0,1). 
Figures Q] (b) and (c) show radiation patterns of the two 
eigenmodes. Figure |TJb) shows the radiation pattern of the 
purely electric eigenmode (p x ,m z ) = (1,0) and Fig. [T]c) 
shows the radiation pattern of the purely magnetic eigen- 
mode (p x ,m z ) = (0,1). Note that both p x and m z ra- 
diate as simple dipoles with a sin 2 9 far field radiation pat- 
tern 5 ^. The two eigenmodes can be selectively excited by im- 
pinging with a plane wave incident along the z-axis with x- 
polarized _E-field (electric eigenmode), or with a plane wave 
incident along the a>axis with y-polarization (z-polarized H- 
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field, magnetic eigenmode). The extinction cross section of 
a single split ring at these two incidence conditions is set by 
cr oxt = 47rfcIm(a £ ; B ) and er cxt = 47rfclm(a_ff ff ). 

Next, we consider extinction and eigenmodes for arbitrary 
values of the cross coupling. It is easy to see that the extinc- 
tion cross section at the two special illumination conditions 
(incident along z, x-polarized, respectively, incident along x, 
with y-polarization) remain equal to o ext = 4:irklm(a.EE) 
and a ex t = 47rfclm(ajjir). However, for nonzero rjc, these 
incidence conditions and polarizabilities do not correspond 
anymore to the eigenvalues and eigenvectors of the polariz- 
ability tensor, which now have mixed magneto-electric char- 
acter. In the extreme case of strongest magneto-electric cou- 
pling (rjc — ^/VeVh), the eigenvectors reduce to (p x , m z ) = 

OMVW^i?) and (Px,m z ) = (1, -iy/m/riE). The asso- 
ciated far-field radiation patterns of these eigenvectors corre- 
spond to coherent superpositions of the radiation pattern of an 
x-oriented electric dipole, and a z-oriented magnetic dipole, 
with a quarter wave phase difference. Figures |TJd,e) show the 
on-resonance radiation pattern, assuming tje — 0.7, r\n = 0.3 
and r\c =0.4 . Note that these parameters are close to the 
limit of strongest possible magneto-electric coupling. Figures 
Q2d,e) reveal that the radiation pattern of each eigenmode is 
non-dipolar. Rather than a sin 2 6 donut-shaped pattern, an 
elongated radiation pattern occurs, with maximum extent in 
the y-direction. The polarization in the far field is linear for 
directions along the cartesian axis, but is generally elliptical. 



B. Extinction cross sections to measure polarizability 

Figure [2] shows the extinction cross section predicted by 
our point scattering model of a single split ring for different 
incidence conditions. In Fig. |2fa), the incident wave vector 
is swept from the z-direction to the y-direction, while main- 
taining x— polarized light. For this set of incidence con- 
ditions the resulting extinction cross sections only depend 
on ccee and ana, and are entirely independent of the off- 
diagonal coupling strength oleh- The cross section increases 
from (Text = 4:TrklmaEE as the split ring is only driven 
by the incident E x field when light is incident along z, to 
Ccxt = 4:irk(ImaEE + Ixnoffjj), as the split ring is driven 
by the incident E x field plus the incident H z field. When the 
wavevector is rotated to the the x-axis, the extinction cross 
section diminishes to AirklmaHH, as the split ring is only 
driven by H z . The chosen values tje — 0.7, t]h = 0.3 
and rjc = 0.4 that we also used for Fig. [Tfd,e) yield ex- 
tinction cross sections cr cx t = 4:irkhna E E = 0.29 /im 2 and 
fcxt = iirklmaHH — 0.13 /im 2 . The predicted er ox t = 
Ankhna.EE = 0.29 /im 2 is consistent with the measurement 
(o"ext = 0.3 /im 2 ) reported by Husnik et al»ii. It is impor- 
tant to note that measurements along cartesian incidence di- 
rections and with linear cartesian polarizations yield only the 
diagonal elements of the polarizability tensor. Indeed, the pro- 
posed measurements form a redundant set of measurements of 
ctEE, ctHH, and (o.ee+(Xhh), but do not provide any insight 
into the magneto-electric cross coupling in the electrodynamic 
polarizability tensor^ 



In order to measure the eigenpolarizabilities, it is neces- 
sary to selectively address the eigenvectors of the polariz- 
ability tensor. As noted above, the eigenvectors in the case 
of strong magneto-electric coupling rjc ~ ^/VeVh tend to 
(Px,m z ) = (l,iy/r) E /rjH) and (l,-i v / r] H /?] E ). These 
eigenvectors require simultaneous driving by E x and H z , with 
a quarter wave phase difference. We note that such fields 
can be generated by circularly polarized light with incident 
wave vector constrained to the xz-plane. Indeed, at maxi- 
mally strong magneto-electric coupling and tje = rjn, cir- 
cularly polarized light incident at 45° from the z-axis would 
selectively excite exactly one eigenmode. Therefore, we ex- 
pect angle-resolved extinction measurements for oppositely 
handed circularly polarized beams to reveal the eigenpolar- 
izabilities. Figure |2jb) plots the extinction cross section for 
right handed circular polarization, as a function of angle of 
incidence in the z-plane, for illumination tuned to the LC reso- 
nance frequency. Naturally, at normal incidence the extinction 
is exactly half the extinction obtained for linear polarization, 
as a consequence of the fact that E y does not interact with 
the split ring at all. Strikingly, the extinction cross section is 
predicted to behave asymmetrically as a function of incidence 
angle. The extinction increases when going to positive angle 
and decreases when going to negative angle. Changing hand- 
edness is equivalent to swapping positive and negative angles. 
A detailed analysis shows that the maximum in extinction cor- 
responds to the largest eigenvalue of the polarizability tensor 
(&ext = 27rklma + ), while the minimum in extinction corre- 
sponds to the smallest eigenvalue {a ext — 27rfclma_). There- 
fore, circularly polarized measurements reveal the eigenvalues 
of the polarizability tensor. Combining such circularly polar- 
ized extinction measurements with the measurements under 
cartesian incidence in Fig.|2|a), therefore allows to extract all 
components of the polarizability tensor. In addition to the 
contrast in extinction, the angle at which the maximum cir- 
cular dichroism occurs is a second, independent measure for 
the magneto-electric coupling strength. The measurements in 
Fig. a) and (b) together hence provide full, even redundant, 
information on rjE,rjH and rjc . 



C. Structural chirality 

The results plotted in Fig. |2{b) show that magneto-electric 
coupling in the 6x6 polarizability tensor directly implies 
structural chirality. It is exhilarating that this interesting phe- 
nomenon first reported by 26 ' 30 for the transmission of arrays 
of scatterers is naturally present in the theory. However, while 
previous analysis of structural chirality focused on transmis- 
sion through periodic arrays, we predict that circular dichro- 
ism already appears in the extinction cross section of a sin- 
gle split ring, with a strength set by how close the magneto- 
electric coupling strength is to its limit y/rjE, r\H- The circular 
dichroism in extinction occurs independently of whether there 
is material loss, as opposed to, e.g., asymmetric transmission 
phenomena through arrays, that are claimed to require dissipa- 
tion^. For maximally magneto-electrically coupled systems, 
the smallest eigenvalue is identically zero, implying that such 
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Incidence angle (°) 

FIG. 2. Extinction cross section cr ex t as a function of the illumination angle and polarization. Blue lines represent a ex t for linearly polarized 
incident illumination, while red lines represent extinction for right handed circularly polarized illumination. Panel (a) shows extinction for 
incidence wave vectors ranging from k z to k y to k x . At normal incidence with k along the z-axis, a ex t is a measure for only oiee as E x is 
the only driving field. Increasing the angle to 90° both polarizations E x and H z excite the dipoles in the SRR, so a ex t is a measure for the 
sum of the terms on the diagonal of the polarizability tensor (oee + oihh)- Changing the angle to couple only the H z of the incident light to 
the SRR gives a ex t that is a measure for purely uhh- Panel (b) a ext as a function of the incident angle in the xz-plane (wave vector ranging 
from —k x to k z to k x ). For right-handed circular polarization minima and maxima in a ex t occur as a function of angle, which are a measure 
for the eigenpolarizabilities a_ and a+, respectively. Both sets of measurements in panel (a) and (b) together provide information on all the 
components of the polarizability tensor, q_bb, &hh, and otEH- 



a scatterer is transparent for one circular polarization, and 
achieves its strongest scattering for the opposite handedness. 
We expect that our 6x6 polarizability tensor can be success- 
fully used to describe all structurally chiral scatterers reported 
today, as well as clusters and periodic arrays thereof. 



V. A COUPLED SYSTEM: SPLIT RING DIMERS 

So far, this manuscript has focused purely on the scattering 
properties of single magneto-electric point scatterers. In the 
remainder of the paper we illustrate that our method can be 
easily used to analyze multiple scattering by magneto-electric 
scattering clusters. In order to calculate the response of a sys- 
tem of coupled magneto-electric dipoles, we generalize the 
general self-consistent equation that describes scattering of 
clusters of electric dipoles p as reviewed in 41 . Assuming a 
system of magneto-electric point scatterers situated at po- 
sitions ri . . . rjy, the response upon illumination by an inci- 
dent field (E[ n (r), H- m (r)) is determined by a set of N self 
consistent equations for the induced dipole moments in each 
scatterer. The dipole moment induced in scatterer n with po- 
larizability tensor a„ is 




(29) 



Using this equation we can attempt to reinterpret recent mea- 
surements that evidence significant coupling in split rings in 
2D arrays, as well as in oligomer a 9 ' 23 - 33 - 34 . Here we focus 
on the extinction of a dimer of split rings in socalled 'stere- 
odimer' configuration, first studied by Liu et all 3 - 3 -. Figure [3] 
shows such a 'stereodimer', consisting of two SRRs in vac- 
uum (V — 200 x 200 x 30 nm 3 , resonant at a wavelength 
around 1500 nm), both parallel to the xy plane, vertically 
stacked with a small height difference of 150 nm. The up- 
per SRR is rotated by a twist angle i/j around the z-axis. On 
the basis of the report by Liu et al. 33 , we expect two reso- 
nance peaks with an angle dependent splitting, which can be 
explained in an LC model as the summed effect of electric 
dipole-dipole coupling and magnetic dipole-dipole coupling. 

We calculate the extinction versus twist angle and wave- 
length of an incident beam incident from the +z direction, 
with x-polarization. This beam directly excites p x in both 
rings, which also drive each other. We first analyze the experi- 
ment assuming that there is no magneto-electric coupling term 
(setting ijc =0, although we keep tje =0.7 and r\n — 0.3). 
As Fig. [3b) shows, the extinction shows a single strong res- 
onance that is blueshifted relative to the single SRR reso- 
nance at 200 THz. As a function of twist angle, this broad 
resonance redshifts to 200 THz at a twist of 90°, and shifts 
back to 220 THz at a twist of 180°. There is no sign of a 
second resonance, which might be hidden below the strong 
resonance. To bring out the second resonance more clearly, 
we reduce the loss in Fig. |3jb), to a 10 times lower value 
7 = 1.25 • 10 13 s _1 ) for gold in Fig. (c) and (d). For this 
almost absorption-free system, Fig. [3J C ) indeed shows two 
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FIG. 3. Extinction cross sections a ax t versus frequency and twist angle for an SRR stereodimer structure. Panel (a) shows the geometry (top 
view and side view) in which two SRRs are vertically stacked. The upper SRR is rotated around the z-axis by the twist angle ip. We calculate 
extinction for light impinging from the ^-direction with polarization along x, i.e., along the base of the lower SRR in (b,c,e,f). In (d,g) we 
use 45° incidence in the :cz-plane, so that the H-field of the excitation light directly couples also to the magnetic polarizability. Panels (b), (c) 
and (d) show extinction assuming no cross coupling term (rjc ~ 0) while (e), (f) and (g) show extinction assuming strong magneto-electric 
coupling r\c = 0.4). Panels (b) and (e) assume the damping rate of gold 7 = 1.25 x 10 14 s _1 . To more clearly bring out the four mode 
structure, we reduce the damping ten-fold for the calculations in (c,d,f, g). There are four modes present in the system. White lines in (b,e) 
indicate the frequencies of the modes, as taken from the resonances in the low-damping case, i.e., the resonances in panels (d,g). Since 



resonances in extinction. The blue shifted resonance is now 
observed to cross with a narrow red shifted resonance. The 
crossing is symmetric around 90° and is consistent with the 
hybridization of an electric dipole fixed along x, with a sec- 
ond one above it twisted by an amount t/j. The two branches 
have a very different width and strength, consistent with the 
fact that a symmetric configuration of dipoles couples more 
strongly to external fields (blue shifted resonance), than an 
antisymmetric 'dark' configuration (red shifted resonance). 

To verify whether the two resonances observed in Fig.[3fa) 
are all resonances in the system, we change the angle of in- 
cidence to 45° in the xz plane, so that the exciting field has 
an H z component to drive the split rings, in addition to an 
E x component. Figure[3jd) shows that in this case four reso- 
nances occur in extinction. In addition to the two curved bands 
excited by E x , there are also two non-dispersive bands with a 
twist independent splitting. Obviously, these bands are due 
to the coupling of two magnetic dipoles in symmetric (broad 
and intense band) and antisymmetric head-to-tail configura- 
tion. The existence of four instead of two modes is a new 
insight compared to LC circuit models 33 ' 36 , but is logical in 
view of the fact that split rings have both a magnetic and an 
electric response, which are decoupled under the assumption 
i] C = 0. 

Next we analyze the extinction in presence of magneto- 
electric coupling, setting ?yc = 0.4. Again, we first examine 
the extinction in presence of realistic loss (7 = 1.25-10 14 s~ x ) 
for gold in Fig.[3]e). As also predicted by FDTD simulations 
by Liu et al. 33 , there appear to be two bands. The blue-shifted 
band is again very broad, but now has a frequency shift away 



from the single SRR resonance that is significantly larger for 
twist angle 180° than for 0°. These effects were explained 
by Liu et al. as due to an additive (subtractive) correction to 
the dominant electric hybridization at twist angle 180° (0°) 
that occurs due to magnetic dipole coupling. A surprise is that 
the diagram is not symmetric anymore around 90° twist as in 
the case of zero magnetic coupling. Instead, the extinction ap- 
pears to show an anticrossing at twist angle 60° These features 
were also predicted by FDTD simulations by Liu et al.— How- 
ever, the presence of an anticrossing at twist angle ip = 60° 
could not be interpreted Liu et al 33 within an LC electrostatic 
circuit model, except by invoking higher order multipolar cor- 
rections. Here we see that a purely dipolar model may also 
explain all features of the experiment provided that magneto- 
electric coupling is accounted for. While we do not claim that 
multipolar effects are not present in actual experiments, it is an 
important insight that split ring polarizabilities with magneto- 
electric coupling terms may provide much richer physics then 
expected from electrostatic circuit theory. A main advantage 
of point dipole theory is that the underlying mode structure 
does not need to be recouped from FDTD simulations, but is 
easily resolved by repeating a calculation of extinction cross 
sections with low loss (as done in Fig. [3]), or by analyzing the 
poles of the coupling matrix in Eq. d29l ) that relates (p, m) 
to (-Em, Hi n ). The computational effort for N split rings is 
equivalent to diagonalizing or inverting a 6A x 6 A matrix. 

To more clearly bring out all the resonances we artificially 
reduce the damping 7 = 1.25 • 10 13 s _1 to ten times less than 
the damping of gold, and plot the response of the system un- 
der normal incidence (f) and 45° incidence (g) in Fig.[3](f,g). 
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The anticrossing at twist angle ip — 60° appears to be due to 
the coupling of four modes, as opposed to the intuition from 
LC circuit theory that only two resonances anticross. The ex- 
istence of four, rather than two modes in a split ring dimer ap- 
pears surprising and is a second indication of the rich physics 
of magneto-electric scatterers. Intuition from LC circuits is 
that although the subspace of driving fields is two dimensional 
(E x and H z ), nonetheless only one mode per split ring exists. 
The usual reasoning in LC models is that the relation between 
electric and magnetic dipole moment is completely fixed and 
independent of driving, since the loop current and accumu- 
lated charge are directly related. Such a constraint is not gen- 
eral: in electrodynamic multipole expansions, magnetic polar- 
izabilties are determined independently from the electric ones. 
The intuition from LC theory that there is only one mode per 
scatterer is only retrieved in our model right at the limit of 
strongest magneto-electric coupling rjc = ^/VeVh, since in 
that case one polarizability is identically zero. We note that 
the values tje — 0.7, ijh = 0.3, rjc = 0.4 used in this work 
(that we fitted to our angle-resolved transmission experiments 
on 200x200 nm Au split rings on glass) are close to the limit 
of strong magneto-electric coupling. Whether a general argu- 
ment exists why physical scatterers are or are not exactly at the 
limit of strongest magneto-electric coupling rjc — Jf]Ef]H is 
a question outside the scope of this paper. 



VI. CONCLUSION 

In conclusion, we have developed a new multiple scattering 
theory by means of which we can calculate scattering and ex- 
tinction for any magneto-electric scatterer with known polar- 
izability tensor, as well as for arbitrary finite clusters. As op- 
posed to LC circuit models, our new model obeys energy con- 
servation, contains all interference effects, and allows quanti- 
tative prediction of absolute cross sections, spectral linewidths 
and lineshapes. While outside the scope of this paper, the the- 
ory is readily extended to deal with arbitrary periodic lattices 
by generalizing Ewald lattice sums 41 to deal with both E and 
H. Since the electrodynamic polarizability tensor can be di- 
rectly constructed from electrostatic circuit theory, we expect 
that our model is readily applicable to many current experi- 
ments using chiral and nonchiral metamaterial building blocks 
for which electrostatic models have been proposed. 

Our model does not give any insight into whether the re- 
sponse of a given structure is truly dipolar or not. Also, our 
model does not provide any insight or quantitative predictions 
based on microscopic considerations for the magnitude of the 
polarizability. For such microscopic considerations, based on, 
e.g., current density distributions derived from full wave sim- 
ulations, we refer to— ~ 20 ' 35 i 37 ' 38 . Rather, our model allows one 
to verify if specific data or microscopic calculations are con- 
sistent at all with point dipole interactions, allowing to ver- 
ify or falsify common intuitive explanations in literature that 
have sofar always been based on electrostatic considerations. 
Also, our model allows one to assess if a single polarizability 
tensor indeed can describe a range of different experiments 



with, e.g., split ring clusters, as should be expected from a 
consistent model. Finally, our model is the simplest electro- 
dynamical model to consistently describe how metamaterials 
and photonic crystals are formed from magneto-electric scat- 
terers. A first step is to confirm the parameters used in this 
work for r\E-,r\H and r\c by targeted experiments. While the 
value for t\e used in this work is consistent with the extinction 
cross section measured by Husnik et aljii, we propose that 
the new insight that magneto-electric coupling is far stronger 
than the magnetic polarizability be confirmed by off-normal 
circularly polarized extinction measurements as proposed in 
section [TV] 

The most important property of our theory is that a polar- 
izability tensor validated for a single scatterer can readily be 
used to predict all quantitative scattering properties of com- 
posite lattices and antennas. We hence expect that new in- 
sights can be obtained in effective medium constants of meta- 
material arrays. Our analytical model not only facilitates de- 
sign, but will also for the first time allow to determine rig- 
orously whether, even in the ideal case (no loss, no multi- 
pole corrections), metamaterial building blocks can give rise 
to a desired e and fj,, despite the large importance of elec- 
trodynamic correction a 7 ' 9 ' 60 . In addition to generating new 
insights for metamaterials, our theory also opens new de- 
sign routes for gratings and antennas with unprecedented po- 
larization properties. As an example, in this paper we an- 
alyzed the four mode anticrossing due to magneto-electric 
coupling in stereo-dimers. This analysis is easily extended 
to magneto-electric Yagi-Uda antennas, diffractive gratings 
of chiral building blocks, and magneto-inductive waveguides 
that may provide new ways to control the propagation and 
emission of ligh t 46 ' 61-62 . 
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Appendix A: Unit system 

Throughout this paper we used units that significantly sim- 
plify notation throughout, as they maximimize the inter- 
changeability of electric and magnetic fields. Conversion to 
SI units is summarized in Table |TJ For the conversion in Ta- 
ble UJ we use e for the host dielectric constant, c for the ve- 
locity of light, and Z for the impedance of the background 
medium. In this unit system, a plane wave has \E\/\ H\ = 1, 
and intensity I = \E\ 2 / (2Z), since the Poynting vector is 
S = l/(2Z)Re(E* x H). In these units, the cycle-averaged 
work done by an electric field E to drive an oscillating p 
equals W = 2irk / ZIm(E ■ p). The magnetic counterpart 
is W = 2-Kk/ZIm(H ■ m) 
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Quantity 


Symbol 


Relation to SI 


Electric field 


E 


Est 


Magnetic field 


H 


ZHsi 


Electric dipole moment 


V 


p S i/(47re) 


Magnetic dipole moment 


m 


msi(Z/(47r)) 


Electric-electric polarizability 


OLEE 


a| s /(47re) 


Magnetic-magnetic polarizability 


OLHH 




Electric-magnetic polarizability 


OLEH 


a| I H (c/(4 7 r)) 


Magnetic-electric polarizability 


OLHE 


aj} E (Z/(4n)) 


Electric-electric Green tensor 


Gee 




Magnetic-magnetic Green tensor 


Ghh 


^G% 1 e 


Electric-electric Green tensor 


Geh 


Att/ZG^e 


Magnetic-magnetic Green tensor 


Ghe 





TABLE I. Conversion between SI units and the unit system used 
throughout this paper. 
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